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FROM TRANSIENT-RESPONSE DATA 
By Marvin Shinbrot 


SUMMARY 


The problem of calculating airplane stability parameters from the 
aircraft response to an arbitrary disturbance is considered. To calcu- 
late the coefficients of the linear differential equation which describes 
the airplane transient response, application is made of a classical 
least- squares curve fitting method. 

It is shown that the method is applicable, although somewhat, 
cumbersome, when the input is an arbitrary function of time. Certain 
inputs are demonstrated to lead to simplification in the application of 
the method. Examples are given illustrating the means of using the 
method and showing its practicability. 

Finally, an appendix, in which Prony* s method and a generalization 
thereof appears, is presented. 


INTRODUCTION 


The determination of the stability and control parameters of a 
dynamical system from its measured response has been a subject of 
increasing importance and interest, both from the standpoint of basic 
aerodynamic research and automatic stabilization of airplanes . This 
problem is essentially one of curve fitting and may be stated as follows: 
Given the time history of an airplane response to a known transient dis- 
turbance, and assuming a certain form of linear differential equation 
with constant coefficients describing the relation between response and 
disturbance, required to find the coefficients of this differential 
equation such that the sum of the squares of the differences between the 
given response and the one corresponding to' the differential equation is 
a minimum. 



2 


MCA TN 2341 


A number of methods for the solution of this problem have been 
advanced by others and are noted in reference 1. However, these methods 
do not apply the least squares principle in the correct sense described 
above . 

The method presented herein involves obtaining the solution of the 
< ^-^ eren 'ki a l equation in a form amenable to rigorous application of the 
method of least squares and such that a completely arbitrary input may 
be treated. A number of means are available for the fitting of a 
function in this form by least squares. The solution by means of a 
Taylor' s series expansion seemingly gives good results and is used 
throughout this paper. 

Considerable simplification of the method is possible when free 
oscillation data are available or when the input can be exactly or at 
least closely represented by a function the analytical form of which is 
known. Solutions for a pulse, a step, or a ramp input are given in 
detail or are indicated. The general solution for an arbitrary input is 
also included. 

■ ' # 

Selected examples of different inputs applied to the same physical 
system (airplane) are presented, which illustrate the different means of 
applying the method. 


The report has been so organized that the engineer who is not 
interested in the derivation of the formulas used may read the first 
section of the report entitled "Statement of the Problem" and proceed 
from there directly to the section on examples. 


METHOD OF ANALYSIS 


Statement of the Problem 


f Consider a quantity which has been measured at a set 

\ fc ip i = 0, 1, . . . , V— 1 of V values of t. Suppose the initial 
conditions on q m , obtained from the experimental data, are 


9m(0) - Qq 



where n is determined below. 
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A physical interpretation of these quantities may he had by con- 
sidering the data q m (t) as a time history of the pitching velocity of 
an airplane in response to an elevator deflection F( t) . 

Consider next the differential equation 

Po(D)q( t) = P X (D)F( t) (2) 

where D is the operator d/dt, P 0 = D 11 + D 11-1 + . . . + a Q and 

Pi = Cm?) 111 + c m-i D 10-1 + • . . + C Q are two polynomials in D, 

and F(t) is a known forcing function. Consider the set of all possible 
solutions of this differential equation obtained by varying the constants 
at and Cj. Let q c (t) be that solution of equation (2) for which 

v-1 

M = £ Cq(ti) - q^ti )f (3) 

i=0 

is a minimum, and subject to the same initial conditions as q m - 



It is then desired to find those values of the a± and the Cj which 
correspond to q c . These values will be unique, since to each set of 
these constants there corresponds one and only one solution of the differ- 
ential equation if the initial conditions are determined. 

In most practical problems , m is less than n, and this assump- 
tion will be made throughout this paper. The extension of the method 
described herein to the case where m is equal to or greater than n 
should be clear. 
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I'et Xi, X 2} . . . , X n be the zeros of the polynomial P 0 (D). It 
is now assumed that f Xj, i / j. The extension of the following 
solution to the case where not all the Xi are distinct can be found 
in any textbook on elementary differential equations and needs not be 
included here (reference >2) . It is well known that if. 


Po'(D) = ^ P 0 (D) 
the solution of equation (2) is 


q(t) 





F(T)dT 


} 


(4) 


where the constants Aj are functions of the initial conditions q(o) , 



and of the constants C x , C 2 , . . . , C m . The method described in this 
paper consists of fitting q(t), by a least squares procedure, to a 
function of the form ( 4) . Several methods for fitting such a function 
may be found in the literature. In this report the classical device of 
linearization and iteration by means of a Taylor’s' series (reference 3 , 
p. 214) is used. Therefore, a detailed study of this linearization 
follows, but it must be understood that it is in no way essential to the 
method. Fitting by steepest descents (reference 4), for example, could 
be used in place of the Taylor’s series iteration. The mechanics of this 
iteration are cumbersome in the. general case. They will, therefore, be 
described for particular cases, thus, it is hoped, making the method 
clear for any special case which may arise. 

Formulation of the Method When Free 
Oscillation Data Are Available 


The method .— Suppose there exists a T such that F(t)=0 for all 
t > T. Then, for t > T, the expression 


.•u.i) L- 


A 1+ ^ 

P ° ^O 


e F(r)dT 


occurring in equation (4) is constant. Call this constant 



MCA TN 23^1 


5 


Then for t > T 




(5) 


By some means (see e . g . , reference l) , first approximations to the 
constants occurring in equation (6) are found. Let B 2 (°), 

. . . , B n ^°^ Xi(°), X 2 (°), . . i , X n (°) he the approximations to 
B 1# B a , . . . , Bn, Xi, X 2 , . . . , X n > respectively. The function 
q(t) is now expanded by means of a Taylor series about the "point" 


* 

Bi(°), . . . , Bn(°), X r (°), 



with all terms of order higher than the first omitted. From equation (5) 


dq 

dBi 


X^t 


dq X^t 

} t ~ ^ e f i = 

o\± 


n 


Denote the values of q, ^L, at the "point" B^ 0 ), . . . , Bh^°^, 

dBi dXi / \ . L / \ 

. , X n (°) by q(°), * respectively. 


,(°) 


Then by Taylor's theorem, omitting all higher order terms. 


Aq 


where 


- x r^V 0) ^ ^v o) 




i=l 


I 


i=l 


\^Xi/ 


AXi 


Aq = q - q(°) 

ABi = Bi - Bj.(°) 
AXi = Xi - X^ 0 ) 


(6) 


i-1, • • • y n 
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Therefore, from equation (6) 


41 = 


n 


£ 6 1 * + Bt(°) t 4\il 

1=1 J 


(7) 


The desired equations which lead to the minimization of 
v-i 

12 


M= Y q c(ti) - q m (ti) 


i=0 


may now be found as follows: Let 

^m = <lm “ 


Then 


v-i 

M= I 


4q(ti) - 4q m (ti) 


i=0 


The equations for the increments then become 


dM 


^(ABi) 


dM 




= 0 


= 0 


i=l, . . . , n 


>( 8 ) 


Since the expression ( 7 ) for Aq is linear in (ABi) and (AXi) 
the equations (8) will be precisely those found when fitting redundant 

a ^m ear fUnction least squares (reference 3 , Sf 20Q-211) 
Equations (8) are solved for the increments ABi and AXi, and the * 

resulting values are added to and , respectively, to 

give the second approximation. 

The iteration process by means of Taylor’s expansions is now 
repeated until two successive iterations give the same values for the 
parameters to the desired number of significant figures 
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It is of course true that in ma n y cases, especially those in which 
very good data are available, no such Taylor series expansion is neces- 
sary, Prony' s method (appendix A) giving sufficiently accurate values of 
the parameters. Examples have been found, however, where Prony' s method 
did not lead to sufficiently good values, and the Taylor series was 
essential. 

Final calculation of the coefficients o f P 0 (D) and Pi(D).- All 
that now remains for the problem to be completely solved is the computa- 
tion of the coefficients of P 0 (D) and Pi(D) from the calc ula ted 
values of Bi and Xi. 

Since P 0 (D) = D n + aj^ D 11-1 + . . . + a 0 ,~ and since Xi, . . ., X n 
are the zeros of this polynomial 


a n— 1 ~ 


a n— 2 - 


a o = (-i) n n Xi 

i=l J 

(reference 5, p. 29). 

The constants Bi were defined as follows: 


n 



i=l 


n 

y *i\j 

/ ( 9 ) 

i< J 


B i = Ai + 


PiUl) 

Po’Ui) 



e XiT F( t ) d T 


where T is that value of t where free oscillation begins. It is 
well known that the Ai may be solved for as functions of the constants 
Ci> . . . , C m . The integrals occurring in the definition of Bi can 
now be found graphically since Xi and F(t) are known. There are, 
therefore, n equations in the m-1 unknowns C Q , . . . , C m (the 
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fore, JieaSr^T 4 !^" £^1 £“» 

r^irr^ r r 

PpMo«?!): e “ w "»«“ -»• ^ 


The Method for a Second-Order System 


a „ £ Ca ® e 0f g f eat Poetical interest in aerodynamics is that of the 
second order, and in the present section this casHm he treated tl 
detail. A second-order system is one in. which equation (l) becomes 


(D e - ai D - a 0 ) q(t) = (C X D - C Q ) F(t) 


An example of such a second-order system which occurs in fiP rn_ 

r 7 be ° btained ^ letting q be the response in pitching 
velocity to an elevator input F(t). in this case (reference l) 8 


ai = i* - & - & 

mV 0 Iy I y 

a 0 =-ia x ^-^ 

mV o iy Iy 


C 0 =^xi2L-^xi 

iy mV o Iy mV 0 

Cl = & - is. X & 

Iy mV 0 Iy 


where 

m mass of airplane 

V 0 trim velocity of airplane 

Iy pitching moment of inertia 
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L lift force 

M pitching moment 

5 elevator input F( t.) 

a angle of attack 


da 

dt 


1*5 


la 


M 8 


hL 

38 

3l 

3a 

3m 

S& 


Ma 


3m 

3a 

3m 

34 


Mq 


3m 

3q 



The general solution of (D 2 + a x D + a 0 )q = 
obtained, from equation (4) which becomes, if X x 
of x2 + a x x + a Q = 0, 


q = A x e Xlt + k 2 e Xst + c ^l* c Q e Xlt 

Xi-X.2 



C X X2+C 0 X,2t 

e 



e X2T F(T)dT 


Differentiating, an expression for 4 
letting t=0, it may be seen that 

A x + A 2 = q(0) 


dq 

dt 


a x \ x + A 2 X 2 = q(o) — 


(C X D + C 0 )F may be 
and X 2 are the roots 


e XlT F(t )dr + 


may be found, and 


C x F(0) 
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or, using equation (la). 


Ai 


\l“"X 2 


a 2 


X 2 — ! ^1 


Therefore, 


q = 


q 0 -q 0 i2^iF 0 + ( c i^i+c 0 ) fo e XlT F(T)dT 


Xi—X 2 




4o-^o^i-CiMCxX g +C 0 ) fl e XgT F(r)dr ^ 

Xa-Xi 


(5a) 


Pulse input .— When a pulse— type input has "been applied to a second- 
order system, that is, when there is a T such that F(t)=0 for all 
t > T, 

' * ; 1 

Bl * Al * / e ~ XlT F< T)iT i 

j> (10) 

- A = ♦ fj °~ X2T F(T)dT 

vising the same notation as before. Then, for t > T, 

q = Bi e^i* + B 2 e x 2 t 


Suppose now, as is generally the case, that Xi and X 2 are con- 
jugate complex numbers. The form previously used for the description of 
the general method may still he followed in this case; however, simplifi- 
cations occur if this notation is abandoned. If and X 2 are con- 

jugate complex, B x and B 2 must be so also, since q is real. Let 
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\ x = Z + Z'i 
Bi = p + p'i 

Then 

X 2 = Z - Z'i, B 2 = p - p'i 

and 


( 11 ) 


q = 2e Zt (p cos Z't - p’ sin Z*t) (12) 

It is assumed that at this point, first approximations 2 0 , Z 0 ’ , 

P 0 , 3o' "^e constants Z, 2', p, p', respectively, have "been 

obtained. (See appendix A for one method for finding such approxima- 
tions.) 


The formulas . needed for the Taylor's series iteration vill now he 
derived. From equation (12), 

s 2te*^ (p cos Z't — p' sin Z't) 

s — 2te^ (p sin Z’t + p* cos Z't) 

dz» 

* ‘ / U3) 
= 2e** cos Z't 


~ = - 2e zt sin Z't 
dp' 


J 


Using a 
the Taylor's 


similar notation to that used in the previous section on 
series expansion, 


v — 1 


m=^ M(ti) - /2a m (ti) 


i=0 


V— 1 


■I i 

i=0 


{ -sr — ) AZ 

\b\J 


<lj 0) 


AP + 




+ q 


(o)_ 


*Zm 


Taking, the derivatives of M with respect to AZ, AZ' , AP, AP' and 
setting them equal to zero will result in the equations, the solution of 
which will give those values of the increments which minimize M. 

These equations thus become 
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As before, it is now assumed that this iteration is repeated until 
two successive iterations give the same values for the parameters to the 
desired number of significant figures. 
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The parameters l, V , p and p* have thus heen determined. The 
solution to our problem will he complete if we can, from these values, 
calculate ao> ai, Co> Ci. From equations (9), 


— Xi + X 2 = 2 1 
a 2 = XiXs = I s + l' z 


A means of finding Ci and C 0 must now he found. In order to 
accomplish this, consider 

* 

rp rp 

J F( t ) cLt = jT e~* T cos Z ! t F(r)dr 


Let 


T 

— i P e sin Z’t F(>r)dT 

J o 

7 = f e~^ r COS l* T r(t)dT 


Let 


Then 


cr* = J' e sin Z't F(x)dT 


y 1 e XlT F(t )dr = a - a*! 


and since \ 2 is the conjugate of \i. 


/ e“ x 2 T F(t )dr = a + c 1 * 


Recalling that B* = p + P’i, it follows from equation (4a) that 
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3 = 


( V ff-Zg* )Ci-cr t Co+Z l q 0 
21 ' 


p» _ HcrZo [ lo+l 1 <J'~F( O) ]C X — oCq 


y (15) 


These equations may "be solved for C Q and c 


1 • 


T he method vhen the forcing f u nction continues throughout the 
-- Otion . It is nov supposed that in the given data, there is no T such 

f ° rC i ne $t ctloa F ( t)=0 for a11 t >T. In this case, the 
coefficient of e J in the expression for q(t) given by equation (4) 

lo nger constant. Prony's method cannot, therefore, he used for 
the first approximation. The simplification which occurs in the Taylor' s 
expansion in the previous case due to the fact that q(t) is a sum of 
exponentials with constant coefficients is also not present here. 

There are cases, however, even when free oscillation does not occur 
vhere the present method can he applied with little difficulty. The 
best examples of such cases are those in which F(t) is known to have a 
certain analytical form. As an example, the two degrees of freedom 
sys em (D + hD + k) q(t) = (CiD + C 0 ) F(t) may be considered, where 

F(t) 


{ 2 * t > 0^ so_ca, TI e< I "step" function. 


For this input, equation (4a) becomes 


q = ^ 1^0 ^-i^23o+Cq gX^t 
^i( ^i~ ^2) 


^2io~ liX.aq 0 +C 0 


lot 


e *2 


Xj_X2 


(4a) 


where, since dq/dt is discontinuous at t=0 if F(t) is a step, 

-t ilgp , dt ’ 0+ indica ting that t is to approach zero through 


positive values only. The apparent anomaly that q 
c i is resolved in appendix B, where it is shown that 


is independent of 


4(0+) = q( 0— ) + C x F(0+) = q( 0 — ) + C x 
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First approximations to the constants 


5 _ Xiq( 0+)— XiXgq 0 +Co _ Xxq( 0-)— XxX2qn+CxXx+Co 

Xi(Xi— X 2 ) Xi(Xi— X 2 ) 

_ X 2 q( 0-*-)— XiXgq 0 +C 0 _ X 2 q( 0-)— XxX g q 0 +C xX g +C Q 
X2 ( X^X 1 ) X g ( X^Xx ) 


B 3 


Co 

XlX 2 


which occur in equation. (4b) can be found by a simple extension to 
Prony‘s method. (See appendix A.) A Taylor^ expansion may then be 
applied to obtain a closer approximation. 

A very important input, closely related to the step is the so-called 
"pseudo— step." In any physical case, an exact 'step can never be obtained 
for an input. The input will always have a certain finite slope near 
t= 0 ; there may also be a certain amount of overshoot or, reciprocally, 
the input may undershoot its steady state value. Any input which rapidly 
(but not instantaneously) attains a constant nonzero value will be called 
a pseudo-step. Suppose then that F(t) is such a function. Then there 
is a value T such that F(t) is a constant for all t > T. Suppose 
this constant is c. Let this function be applied to the second— order 
system (D 2 + ax + Sq) q(t) = (CxD + C 0 ) F(t). The response in q is 
then given by equation (4a). It is here assumed for simplicity that the 
initial conditions q(0) and q(0) are both zero; F(o) is also zero. 
Then 


q(t) 


CxXx+Cp 
Xx— X 2 



e _XlT F( t ) d t + 

Xs>— Xx 



e~ XsT F(r)dT 


For t > T, 


q(t) 


CxXx+Co Xxt 

e 



0 “XlT 


F( T ) d T 



F(t )d t 


+ 


CiX 2 +Cq e \ 2 t 



e -X 2 T F(T)d T + 



e~^2T F( T )dT 
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e XiT P( r)dT 
e _x l T F ( T ) dT 


Is constant and since 
may be evaluated. 


F(t) 


is constant for t > T, and 


_c_ CxXi+Cq _ n c CxX 2 +C 0 
Xi X1-X2 Ux ' X 2 \Pkl ~ Gz > 


CxXi+Cq 

Xx-X 2 J 


e x i T F(x)dr 


Hi, 


C1X2+C0 



e -X 2 T F ( T ) dT _ Ha 


Then Gx, G 2 , Hi, H 2 are constants, and 


q(t) = Hi e x i t + Gx £e x i( ’ b-T )— ij + h 2 e x 2 ^ + q 2 ^e x s( 1 


Thus q(t) is the sum of two exponentials plus a c ons tant , First 
approximations may now be found by the extension to Prony’ s method given 
in appendix A, and Taylor’s series expansions may be used to improve 
these values. 


Another simple example which may be considered is that where 

F ^L = Lt! t>0* In this case the response of the second-order system 
would be 1- ’ — 


<1 = 7ie Xlt + 7 2 e Xst - (nXi+7 2 X 2 ) t - (7x+7 2 ) 


if 


q(0) = q( 0) = 0 


vhere 


Cx^x+Cq _ CiXs+Cq 
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A first approx ima tion may be found by several means, the derivative 
method (see example 3) being only one of these. Again, Taylor's series 
may be applied to improve these values. 

In the most general case, where F(t) is only known graphically or 
tabularly, the method can still be applied; there is a certain amount 
of added labor involved, but in case such an example should occur, the 
method for this general case is outlined below. It Is to be noted that 
this general case is better suited to computation with a high-speed 
electronic computer than it is to manual computation. However, it may 
certainly be applied manually. 

The entire method proceeds from equation ( 4 ). It is assumed that 
the given differential equation has been solved for q in this form. 

Step 1 .— By some means', a first approximation to the parameters 
must be found (several such methods may be found in reference 1) . 

Step 2 .— The function q(t) is expanded in a Taylor's series, all 
teams of order higher than the first being omitted. Herein lies most of 

the computation, since f Q e 1 F(r)dT must be found graphically as a 
function of t. 

As an example, consider again the second— order case (D 2 + axD + Sq) 
q(t) = (CiD + C 0 ) F(t), where F(t) is given graphically.. It is 
assumed for simplicity that 


F 0 = q(0) = q(0) = 0 


Then, using the same notation as before. Ax = A 2 = 0, and q(t) . is 
simply the particular integral: 


e-^ T F( T )d T + ^e^/ e-^F( T )d, 


The computations (especially when it comes to taking derivatives for use 
in the Taylor series) are greatly simplified by the substitution 
C1I1+C0 CxX 2 +C 0 

Bx = B 2 = — . The problem is then restated as follows: 

\x — I2 Xgr- Xx 

To find the best fit for q(t) to a function of the form 


q(t) = Bx e^f e“ XlT F(t ) dT + B 2 e*** J * e“ X 2 T F(r)dT 
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Minimization of M = £i(q c -q m ) with respect to 3 , 3 ’, 2, V , (where 
3 + 3*i = B x , 3 - 3*i. = B 2 , 2 + 2'i = X x , 2 — 2*i = X 2 ) is equivalent 
to minimization with respect to C x , C 0 , a x , and ao, which is the 
desired minimization. The method now proceeds as before, with the first 
approximations being found and with the application of a Taylor* s series . 

Another means which might be employed in this general case would be 
to first fit F(t) to some suitable function, and then find q(t) as 
an explicit function of t as was done above for the step and the ramp. 


EXAMPLES 


The engineer who is interested primarily in applications may read 
these examples immediately after reading the first section of the report 
entitled "Statement of the Problem." 


Example I - Pulse Input 

Consider the input shown in figure 1(a) applied to the secondr-order 
system 


(D 2 + a x D + a Q )q(t) = (C X D + C 0 )F(t) ( 1 6 ) 


The response is given in figure 1(b). The methods and formulas described 
previously can be used directly for this example. If t > 0.4, 


q = B x e Xlt + B 2 e Xst 


(5a) 


^ 3 t 3 * i, X x — 2 + 2 * i, and t ^ 0 . 4 , 


q 

2 


2t 


(3 cos 2 *t - 3 * sin 2 * t) 


( 12 ) 
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It Is now assumed that by some means first approximation l 0 , 2 0 ’ , 
P 0 , Po* > to If V , p, PS respectively, have been found. 1 Prony' s method 
(see appendix A) gives 2 0 = — 0.9S 2 0 * = 7*02, p 0 = O.368I, 
p 0 » = - 2.6775. 

From equation (12) 



t e*^ (p cos 2’t — p’ sin I’t) 


— t e^ (p sin I'.t + P' cos 2’t) 


e^ cos 2't 


- e U sin 2’t 







Prony* s method, for example, may be used here. Prony* s method is 
actually applied to this example in appendix A. An ordinary period 
and damping analysis (described below), which is shorter than Prony’ s 
method, may also be used. The afore— mentioned period and damping anal- 
ysis proceeds as follows: The minima of q lie on a curve q = K e^, 

and q 2 be two successive minima of q. Then 

q (which may be found from figure 1 (b)). 


K = constant. Let 



As for 2' , 2’ 


is equal to ^ 


Then 


q = 2 e lt (p cos 2’t - p* sin 2’t) 


Fitting at any two points and t 2 

(cos 2'tx) p - (sin 2'ti) P' = - e“ Ztl q(t x ) 

2 

(cos 2 * t 2 ) p - (sin 2 ’ t 2 ) p’ = | e " Zt2 q(t 2 ) 

which may be solved for p. and p’ . This analysis gives 2 3-0.9, 

2* 3 7.0. Letting ti = 1.0, t 2 = 2.0, it may be found that p 3 0.+, 
P’ 32.6. 
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Therefore, by Taylor’s theorem, omitting all higher order terms, 
t(p cos 2’t - p' sin 2’ t) A2 - t (p sin 2’t + p' cos 2»t) A2‘ + 

(cos 2't) AP - (sin 2’t) Ap’ = ^ e~ U - (p cos 2't - p* sin 2’t) 

Referring to table I where circled numbers refer to columns, 

@ A2 - @ A2’ + @AP - @ Ap' = © 

Minimizing M leads to the equations (see equations (14)) 

A2 E (gf -A2’ E @ @ + AP E © ©^P’ Z @ © = £ © © 

-A2 Z ©@+A2’Z @ 2 - AP E0@+.AP’ £ @ @ =- E @ @ . 
A2 Z0'©-A2’E 0 @ + APE 0 2 -Ap’ Z@@ =E0@ 

-A2 E© © + A2’E © APE 0@ + AP’ E © 2 =- £© © 

These equations give 

A2 = — 0.0100 A2' = 0.0052 

AP = - 0.0118 AP' = - 0.0320 
or 

2 = - 0.9200 2’ = 7.0252 
P = 0.3563 P' = - 2.7095 

A second iteration was now applied, and it was found that the increments 
were zero to four decimals. 
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Finally, for computing a x , a 0 , Ci and C 0 

a x = - 2 l = 1.84 

a Q = I s + l' 2 = 50.1998 

The functions e — ^ F(t) cos Z't and e~ ^^F(t) sin I’t are now tabu- 
lated and plotted (table II and fig. 2). These two functions must now 
be integrated to find a and cr* . It is to be noted that in this 
example F(t) is exactly equal to t for 0 < t < 0.2, and to 
(0.4 — t) for 0.2 < t < 0.4, and so the integration may be done 
analytically. The integration may be done graphically by means' of a 
planimeter, however, when F(t) is not so clearly a perfect triangular 
pulse, as is usually the case. A planimeter was actually used in this 
example to give a = 0.00493 and a* = 0.0405. Equations (15) are now 
set up: 

0.280 Ci + 0.00493 C 0 = 38.0690 
0.07189 Ci - 0.0 405 c 0 = 5.0061 

which give Ci =' 134.0 and C 0 = 114.2, to four significant figures. 

The true values of the parameters in the preceding example were 

ai = 1.84 
a 0 = 50.2 
Ci = 134.0 
C 0 = 114.4 

Example II — Step Input 

Consider a "step" input F(t) * t > 0 a PP lled to the s y stem 

described by equation (l6). The response is given by equation (4b). 

From the response shown in figure 3 it may be seen that q(0) = 0, 
q(0-) = 0. Then equation (4b) becomes 
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q(t) 


CiXi+C 0 ^ t CiXa+Co Xpt C 0 
Xi(Xi— X2) X2(X2— Xi) X1X2 


Let 


Ci^i+Co 
Xi(Xi~ X2) 


B = C1X2+C0 

Xs(X2-Xi) 


B3 


= _£o_ 
XlX2 


Then, using the same notation as before, with X x = 2 + 2' i, XI = 2 — Z'i, 
Xi + X2 = &i, X1X2 - ao. 


Co Cjao+CoZ 


2a 0 22* a Q 


B 2 = — + Ci a o+Co2 i 

2&o 2 Z 1 a 0 


Let B a - (S + P'l. Then f> - - S’ = - . But B 3 - ^2 . 

o a ° 22 a o a Q . 

Therefore, B 3 = - 2(3 . Then | = e 2t (p cos 2't - p* sin 2't) - p . 

„ J 1 ”* approximations may now be found by the extension to Prony's 
method described in appendix A or hy a generalization of the ordinary 
period and damping analysis. For simplicity, this period and damping . 
method was used to give 2 0 = - 0.91, 2 0 ' = 7.0, p 0 = - 1.11, p Q « = - 9.7 . 
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The Taylor's expansion is now applied 


= e^ (p cos Z't — p' sin Z't) — p 


t e^ (p cos Z't — P' sin Z't) 


MV- 


it 


— t e (P sin Z't + p' cos Z't) 


5(i)"- eW cos 2,t_1 


ap 


gin ,, t 


Therefore, 


e lt (p cos Z't— p* sin Z't)— pj + j^t e lt (p cos Z't — p' sin Z ' t ) J AZ - 
t e lt (p sin Z't + P' cos Z't)jAZ* + £e ?t cos Z't>-1 AP — sin Z'tjAP' 


Or, referring to table III, (l£)AZ + (l9)AZ' + (l5)AP — (lCyAP' = (21) . 
Again minimizing M, it may he found tnat 


AZ = - 0.0100 
AZ' = 0.0242 


Ap = - 0.0292 
AP' = 0.3110 


which in turn give 
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2 = - 0.9200 
2' = 7.0242 
P = - 1.1392 
P' = - 9.3890 

A second Taylor's expansion was then applied (see table IV), giving 

2 = - 0.9201 
v = 7.0258 
- 0 = - 1.1400 

0’ = - 9.3875 

Another iteration would lead to increments which are zero to four 
decimals. 

For the final calculation of the coefficients, 

= - (Xi+X 2 ) = - 22 = 1.8402 
a 0 = X x X 2 = 2 2 2' 2 =50.2084 
C 0 = 2k0 = 114.4715 

Cl = - ? kzl P ,+c Q* = 134.0072 

k 


It is to be noted that a second iteration such as was applied above 
was hardly important, since the results of the first Taylor's series 
give 
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a x = 1.84 


a o = 50.19 
C Q = 114.34 
C x = 134.00 


the true values being, as before. 


a x = 1.84 


a Q =50.2 


C Q = 114.4 


C x = 134.0 


Example III - Input Requiring the General Method 


Consider finally the input shown in figure 4(a) applied to the 
system 


(D 2 + a x D + a 0 )q(t) = (C X D + C 0 )F(t) - (16) 

Such an input might occur, for example, in a stabilized airplane where 
the pitching velocity q is fed back to the elevator to change the 
input . 

First approximations by the derivative method .- The so-called 
"derivative method" will here be applied in order to find first approxi- 
mations to the desired constants. From figure 4, F and q are found 
graphically as functions of t. Then q is plotted (fig. 5), and from 
this, q is found and tabulated. (See table V.) Rewriting equation (16), 


(q)ai t (q)ao - (F)C X - (F)C 0 = - (q) 
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and referring to table V, 

(g)ai + ©ao - © Cl - ®C 0 = - © 

A least-squares analysis is now applied to this equation, giving 

a x = 1.84 
^ = 50.19 

Ci = 133.89 

C 0 = 114.91 


The high accuracy of these first approximations is due, of course, to the 
excellence of the data in the example. "The leading objection to the 
derivative method obviously is the necessity of finding the derivatives 
graphically, which in many cases leads to errors so gross as to make the 
values of the parameters found in this manner entirely valueless. 


The Taylor’s series iteration .— Since F(t) is continuous at t=0, 
the constants q(0), 4(0), F(0) of equation (4a) may be found as usual, 
by inspection of figures 4 and 5. Since q(0) = q(0) = F(0). = 0, from 
equation (4a), 


q = CjXjjHCo e \it r* e XiT p ( T ) cLt + £l^2+2o e X 2 t P* e *aT P ( T ) dT 

^■1 ^2 J 0 \ z -\ 1 J o 


or, letting Bl = , B 2 = , 

A, X*”^2 A, 2 ^^ x 


, = Bl r e X * T 


f e^ lT F(-r)dr + B 2 e^ 2 ^ e^ 2T F(t)dn 


Using the same notation as before, with \ x = 7 + I'i, Bj. = 3 + p'i, 
f Q e — ^ lT F(-r)dT = a — o' i, it follows that = l — 1*1, 

B 2 = p - p'i, /* e~^ 2T F(r)dr = <7 + cr'i, and that, 
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q = 2 e lt [(3cr + 3'p') cos Z't + (3cr’ — 3' cr) sin Z't] (17) 


It will also be noted that 3 = — ^ > 3' = — — ^ . Letting zero sub- 

2 2 Z 

scripts denote the values of the indicated parameters at the first 
approximations , 

Z 0 = - 0.92 
I’o = 7.02 


3 0 = 66.945 


From equation (17) 
dq 

-^ = 2 e 

hi 


It 


3' 0 = - 0.59 


[ (3cr+3* cr* ) t — (3p+3*p' )] cos Z’t + 




[ (3<t'— 3' cr)t + (3 'p-3p')] sin Z't 


} 


6q 

hi' 


= 2 e lt -\ [(3ff , -3’cr)t + (3’p-3p')] cos Z’t - 


{ 


[ (3cr+3'cr' )t - (3p+3’p')] sin Z’tj- 

— = 2 e lt (cr cos Z’t + cr' sin Z't) 

c>3 

■r-^- = 2 e^ (cr' cos Z't — a sin Z't) 
d3' 


s 


/ (18) 


where 


r 


p = 


t F(t) e cos Z't dt 


p' = J' t F(t) e ^ sin Z't dt 
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The function F(t) is now multiplied by the four quantities 

e lot cos 2' Q t, e“* ot sin Z' 0 t, t e~ lot cos Z' 0 t, t e""* ot sin 2' 0 t, 

point by point, and the resulting products are plotted against t 
(fig. 6). These four curves are integrated as functions of t, giving, 
respectively, cr, cr* , p, and p' (table VI). 


The quantities 
/dq\ Z 0 t 

Wo 6 


<Jo 


(8q/82) Q (8q/dZ’) 0 (dq/d3) 0 

e 2 0 t ' 7+ ' 7 * > 


e z ot 


e 2 0 t 


and 


are then computed from equation ( 18 ). 


Then, using a Taylor’s series and table VI, 


© A2 + © A2' + @A3 + @A3' = @ 
Minimizing M with respect to A2, AZ', A3, A3' gives 

A2 = — 0.00 


AZ' = 0.01 
A3 = 0.17 


A3* = - 0.05 


Due to the inherent inaccuracy of a planimeter, only two decimal places 
were preserved here. 

Thus 


z = - 0.92 
v = 7.03 


3 = 67.03 


3' = - 0.64 
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The desired parameters were finally computed to be: 

a x = 1.84 
a Q = 50.28 

Ci -= 134.06 
C Q = 114.69 

1 

CONCLUDING REMARKS 


A method has been described by which the coefficients of the differ- 
ential equation 


5”a + » + 

at n ^ dt“- 1 ' 


. + ax 3a + Soq 

dt 


an* ap-ij? 

m dt m n>_1 dt®" 1 * 


. + Ci — + CqF 
dt 


can be calculated from the graphical knowledge of q(t) and F(t). It 
' is noticed that the method may become somewhat cumbersome if the input 
F(t) is not of certain types. The input which allows the method to be 
applied most easily is one which goes to zero quickly, giving free oscil- 
lation data. A simplification also occurs when the input F(t) is 
known accurately to be a function of a certain type, such as a step, a 
ramp F(t) = t or, perhaps, something of the form S Aj_e^ 

The method has as its primary advantages first, the fact that the 
correct quantity (the sum of the squares of the differences between the 
calculated and the measured quantities) is minimized. The ordinary 
simple equations of least squares may be used in this minimization, since 
the equations of condition (the redundant or inconsistent equations to 
be solved by least squares) satisfy the restrictions .under which the 
least squares solution is normally derived. That is, the equations are 
linear, and only the right-hand sides are subject to error. The so- 
called "derivative method" used in example III minimizes 
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• • • ^ 

2(q + b<J + kq — CjF — CqF) , a quantity with little, if any meaning. 
Other methods, Prony's, the Laplace and Fourier methods, are subject to 
the same objection. 

Another advantage of the method is that the data are analyzed 
directly in the time plane, thus eliminating the possibility of the 
introduction of errors due to graphical integration (in the case of 
Fourier and Laplace transform methods) or differentiation (derivative 
method) . Where the inputs give free oscillation data there is a certain 
saving of time due to this directness, no prel imin ary steps having to be 
taken before the method can be applied. 

For the usual type of problem which is found, those with the pulse— 
or step-type inputs, the method is at least as rapid as any other which 
gives comparable results. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif., Jan. 17 , 1951. 
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APPENDIX A 
PROXY’S METHOD 


A problem which occurs quite often is that of fitting a function 
q(t) to a sum of exponentials. One method used to fit a function of 
this form is known as Prony’s method and is described below. (See also 
references 1 and 3.) 


It is first assumed that the data are given in tabular form with 
equal intervals of the argument t. Suppose measurements of q(t) are 
aken at instants t 0 , t x , . . . , ty-i, v > 2n, where n is the number 
of exponentials to which q is to be fitted. Then it is assumed that 
k o + k(At), At - constant, k = 0, 1, . . . , v— 1. By a proper shift 
of the time axis, t 0 may be taken equal to zero. Let q k = q(t k ) = 

q(kAt). Fundamental to Prony's method is the following theorem: If 

n 




Bi e 


Xit k 


(19) 


i=l 


k = 0, 
tion 


1 , 


v 1, then q satisfies the linear difference equa— 


qk +n + Qtt-i Vi + q k+n-2 + 


+ Q 0 q k = 0 


( 20 ) 


where Qq, 
tion 


^n— 1 are constants such that the roots of the equa— 


x11 + Vi 


Qo = 0 


Xi(At) 

are e , i - 1, . . . , n . The proof is as follows: 


( 21 ) 


^k+n + Sn-i ^k+n^-i + • • • + O^k = B i e Xi ^ k+n ^ At + 

i=l 


n 



i=l 


e Xi(k+n-i)At + 


V Bl e^ At 

1=1 




e \±kAt 



+ 


+ Qq 
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Certainly there 

XiAt X.oAt 
e > e } • 
these constants 
and, therefore, 
the proof. 

Prony’s method consists of writing down the equations 

^n + ^n— 1 *!n-i + • • • + ^o^o = ® 

1n+i + ^n— l *3n + • • • + Qo^i = 0 

q n+2 + ( ^n — i ^n+i + • • . + Qo<1 2 = 0 


exist n constants Qq, . . . , Q n _ 1 , such that 

. . , e are the roots of equation (21). Choosing 

in this way then makes the bracketed expression above 
*lk+n + ^n-i ^k+n-i + • • • + Qq van ish. This completes 


and solving them by least squares 2 for Qf, i = 0, . . . , n-1. From 
the normal equations obtained from the least — squares process, 

Qo> • • • > Qn-i can found. Then the roots of the equation 

x n + Qn-i x 11-1 + . . . + Q 0 = 0 


are calculated, giving e^i(At) therefore, \ i# 

2 Herein lies one of the objections to Prony' s method, since 


v-n 


7 ^k+n + ^n-i 'Ik+n-i + • • • + 


Vk)' 


k=0 


rather than 


T U-f »i.^) 

kK) i=l ^ 


is minimized, the correct minimization procedure leading, in this 
case, to a forbidding amount of calculation. 
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The coefficient, e^** of Bi in equat ion (19) can "be tabulated 
for each i, since Xi is known. Therefore, can be found by a 

second least-squares procedure. This completes Prony’ s method. 

The objection described in the preceding footnote can be overcome * 
by considering the X^ and B^ found by Prony' s method not as the best 
possible values of these 'parameters, but only as a first approximation 
thereof. A better approximation can then be found by means of a Taylor’s 
series. This method is described in the body of this report. 


Extension to Prony' s Method 

Suppose now that it is required to find some way to fit a function 
to a sum of exponentials plus an (unknown) constant. The solution will 
be presented herein for the case in which the exponentials are two in 
number, but the generalization to a greater number of unknowns will be 
evident . 

To be specific, a quantity q must be fitted to a function of the 

form 

q = B x e x i t + B 2 e*®* + B3 (22) 


q is first fitted to a function of the form 


q = Bi e Xlt + B 2 e Xat + B 3 e X3t (23) 

% 

and the condition that X 3 = 0 is put in later.. Prony' s method is 
applied first to equation (2) . As before q satisfies the difference 
equation 


9k+3 + ^2 ^k +2 + ^1 ^k+i + ^o^k 0 (24) 1 

where the constants Qj_, i=0, 1, 2, 3, sire such that the equation 

x 3 + Q 2 x 2 + x + Q 0 = 0 (25) 

has e x i(^) as roots. But X3 = 0. Therefore, equation (4) has 
unity as a root, and 1 + Q 2 +' Q x + Qo = 0. Eliminating Q q (say) 
between this and equation (24) 
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(^k+2 ~ 9}j) ^2 + (^k+i “ 9k) = (*3k — ^k+3^ (26) 

is obtained. This- equation is now solved by least squares for Q x 
and Q e . The exponentials and are now found to be 

the roots of x 2 + (Q 2 + 1) x + (C^ + Q 2 + l) =0. Prony*s method now 
proceeds as before. 

An evident extension of this method may be used to simplify the 
problem of fitting a sum of exponentials when one or more of the expo- 
nents are already known. 

, Example of Prony 1 s Method 

Suppose the data found in column (3) , table VII is to be fitted to 
two exponentials (see example I of the body of this report). In this 
case, equations (20) become 

© + ©Qi + @Qo = 0 

Solving by least squares, 

"Q 1 E© 2 + QoS© ® = - 2 © © 

H 

.Qi£®©+Qo2© 2 =-£® © 

which give 

Q 0 = 0.8320 
Q x = - 1.3922 

Equation (21) thus becomes 

x 2 - 1.3922 x + 0.8320 =0 


or 
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x = e*- 1 ^) = O.696I + O.5895 i 

Therefore, 

Xi = — In (O.696I + 0.5895 1) 

At 

='10 I" 2n J (0.696l) 2 +(0.5895) 2 + i arc tan - 1 

L 0.6961 J 

= - 0.91 + 7.02 i 


X 2 = - O.91 - 7.02 i 


Since \ x and. X 2 are complex conjugate, so are B x and B 2 . Let 

Bi = p + p'i. Then q = 2 e^ (p cos Z't — P' sin I't). Referring to 
table VII, this becomes (7) 0 - (§) P' = (§> . Solving as before by- 
least squares. 


P = O.368I 


p* = - 2.6775 
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APPENDIX B 
DISCONTINUOUS INPUTS 
Consider the differential equation 


+ « £ 1 a 

dt n ^-1 dt 15-1 


+ ai M + “o* - 




dt 1 


.m 


(2) 


A fundamental question which often arises is the following. Suppose 
F(t) or one of its derivatives is discontinuous at a point. At this 
point all higher derivatives fail to exist. What is the meaning, if anv 
of the differential equation (2) at this point? V ’ 


An instance of this problem occurs in example II of the present 
report. There, F(t) is a step: F(t)~fJ'J>°. At t = 0 

dF/dt and all higher derivatives of F do not exist. 

Throughout the following discussion, let x(0+) = lim x(t) 

v /.x t— >0+ 

where x(t) is any function, and t — > 0+ indicate that t is to 

approach zero through positive values only. Consider again the differ- 
ential equation 

q + ax q + a 0 q = C X F + CqF (16) 


where dots denote differentiation with respect to t. Integrate equa- 
tion (16) between the limits — e and t (e>0) to obtain 


[q(t) - q (— e ) ] + ai[q(t) - q(-e)] 




q(t)dt 


■ Ci [F(t)-F(— e)] 


+ Cq 



F(t)dt 


It is now assumed that /_ e q(t)dt is continuous. (In all examples used 
in this report, q(t) represents the pitching velocity of the airplane, 
making j*_^q(t)dt the angle of pitch which, from physical considered 
tionsj may be seen to be continuous.) 
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J* F(t)at = f F(t)dt = t, 


since 



t < 0 
t > 0 


Letting e — •> 0, 


[q(t)-q(0-)] + a x [q(t)-q(0-) ] + k[e(t)-0(O-) ] = C x F(t) + C Q t 


where 



q(t)dt 


Integrating once more with respect to t, this time from — e to +e. 


[q(e)-q(-e)-2eq(0-) ] + a x [0(e)-0(-e)-2eq(O-) ] + 



0 ( t ) dt — 2e0(O— ) 


-Ci€ 


Again letting € — >0, it is seen that q(0+) — q(O-) = 0, since 
0(t) was assumed continuous. Thus q(t) is continuous at zero. Going 
back to equation (19) and letting t — >0, q(0+) — q(O-) = C x F(0+) = C x . 
Thus, q(t) is discontinuous at zero, and the difference between the 
right-hand and left-hand limits of q(t) as t approaches zero is C x . 
Thus, writing equation (4a) in terms of q(O-) rather than q(0+), 

/j _\ X x q(0— )— X x X2q(0)+C x X. x +C o x, x t X2q(0-)— A x X 2 q(0)+C x \ 2 +C o x 2 t 

q(t) e + e + 

X1X2 


and q(t) may be seen after all to be dependent on C x . 
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TABLE II.- FURTHER CALCULATIONS NECESSARY FOR EXAMPLE I 


3 

4 

5 

6 

7 

8 

9 

10 

0.92t 

e-.92t 

7.0252t 

coa (5) 

am© 

F(t) 

© 

® 

© 

Lo:©:©: 

0 

1.0000 

0 

1.0000 

0 

0 

0 

0 

.046 

.9550 

.3513 

.9389 

.3441 

.05 

.04483 

.0164 

.092 

.9121 

.7025 

.7632 

.6461 

.1 

.06961 

.0589 

.138 

.8711 

1.0538 

.4943 

.8693 

.15 

.06469 

.1136 

.184 

.8319 

1.4050 

.1650 

.9863 

.2 

.02745 

.1641 

.230 

.7945 

1.7563 

-.1844 

.9828 

.15 

-.02198 

.1171 

.276 

.7588 

2.1075 

-.5113 

.8594 

.1 

-.03880 

.0652 

.322 

^7247 

2.4588 

- 7758 

.6310 

.05 

-.02811 

.0229 

.368 

.6921 

2.8100 

-.9455 

.3256 

0 

0 

0 





























TABLE III.- TAYLOR'S SERIES ITERATION APPLIED TO EXAMPLE II 



NACA 



IV.- SECOND ITERATION FOR EXAMPLE II 
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TABLE V.- TABULATION OF F(t), F(t), q(t), <i(t), AND q.(t ) FOR EXAMPLE m 


0 . 

.1 

.2 

.3 

.4 

• 5 
.6 

• 7 
.8 

.9 

1.0 

1.1 

1.2 

1.3 

1.4 

1.5 

1.6 

1.7 

1.8 

1.9 
2.0 
2.1 
2.2 

2.3 

2.4 

2.5 

2.6 

2.7 

2.8 

2.9 
3.0 


• F 

0 

.5893 

.8206 

.6521 

.2251 

-.2290 

5062 

5141 

-.2947 

.0175 

.2694 

.3606 

.2779 

.0869 

1103 

-.2258- 

-.2226 

1221 

.0152 

.1228 

.1582 

.1182 

.0329 

-.0525 

1005 

-.0962 

-.0503 

.0100 

.0558 

.0693 

.0501 


-3-33 
-4. 81 
-3.92 
-1.46 
I.23 
2.92 
3-05 
1.8l 

-.01 

-1.53 

-2.11 

-1.67 

-.57 

.60 

1.31 

1.32 
.75 

-.05 

-.70 

-.93 

-.71 

-.22 

.29 

.58 
.57 
• 31 

-.04 


g 

0 

3.940 
10.939 
12.974 
5.847 
-7.837 
-20.407 
-24.019 
-15.576 
1.300 
18. 214 
26.577 
22.181 
7.319 
-10.555 
-22.775 
-23.785 
-13.711 
2.053 
15. 764 
21.207 
16.504 
4.663 
-8.250 
-16.219 
-16.060 
-8.637 

1.945 

10.555 

13.468 

9.997 


4 

0 

68.4 
57.1 

-23.9 

-113.9 

-146.1 

-91.3 

24.6 

137-9 

184.8 

-138.3 

21.5 
-105.4 
-178.7 
-163.6 

-71.0 

50.6 
l4l.4 
160.2 
103.1 

2.5 

-91.2 

-134.9 

-112.8 

-4l.0 

42.6 

98.6 
104.0 

61.6 
-4.3 

-61.1 


% 

9 5T 
327 
-532 
-978 
-702 
110 
938 
1266 
886 
5 

-895 

-1334 

-1092 

-312 

594 

1171 

1155 

589 

-219 

-866 

-1058 

-74 1 

-107 
520 
849 
755 
- 324 
-211 

-595 

-666 

-426 


z© 2 = 329048.34 z© 2 = 6660.992256 z@© =-0.428921 

z @© = 190.3603 Z@© = -82.75361 £©© = 16887.46 

z ©© = ^691.967 2©@ = +98.2381628 z 0 2 = 2.66225444 

z ©© = 83.67703 z ©© = -334446.428 z©© = -4835.7037 

£©© = 23599.0 z © 2 = +159.8717 ^ ^ 
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TABLE VI.- TAYLOR'S SERIES ITERATION APPLIED TO EXAMPLE III 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

Rov 

t 

a 

0* 

p 

P* 

<lo 

(sfl 

m 

© 

(&)„ 


©-© 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

2 

.1 

.0297 

.0149 

.00191 

.0029 

3.9370 

-.0110 

-.2114 

.02947 

.00711 

3.940 

.0030 

3 

.2 

.0692 

.0884 

.00755 

.0138 

10.9405 

.5410 

-.6242 

.08210 

.04459 

10.939 

-.0015 

4 

.3 

.0526 

.1813 

.00287 

.0309 

12.9966 

1.3609 

-2.35*15 

.09799 

.10445 

12.974 

-.0226 

5 

.4’ 

.0075 

.222 

-.01278 

.0449 

5.8985 

-.0912 

-4.3174 

.04535 

.14683 

5.847 

-.0515 

6 

• 5 

.0094 

.226 

-.01141 

.0572 

-7.7700 

-3.0022 

-3.8860 

-.05688 

.13094 

-7.837 

-.0670 

7 

.6 

.0550 

.270 

.01376 

.0850 

-20.3630 

-5.9425 

-1.5001 

-.15168 

.04685 

-20.407 

-.044o 

8 

.7 

.0682 

.364 

.0218 

.1323 

-24.0440 

-8.0511 

3.5835 

-.18022 

-.07340 

-24.019 

.0250 

9 

.8 

.0270 

.435 

-.0094 

.1757 

-15.6851 

-5.1771 

9.9342 

-.11866 

-.17151 

-15.576 

.1091 

io 

• 9 

.0001 

.447 

-.0317 

.201 

1.1454 

2.3728 

11.6801 

.00683 

-.19527 

1.300 

.1546 

n 

1.0 

.0325 

.465 

-.0006 

.235 

18.0924 

9.6213 

7.8174 

.13400 

-.12840 

18. 214 

.1216 

12 

1.1 

.0700 

.541 

.0385 

.303 

26.5549 

14.3646 

-.1106 

.19833 

-.00064 

26.577 

.0221 

13 

1.2 

.0500 

.633 

.0150 

.383 

22.2894 

12.8773 

-10.8027 

.16760 

.12714 

22.181 

-.1084 

14 

1.3 

.0059 

.670 

-.0400 

.436 

7.5236 

3.1852 

-17.4231 

.05790 

.19405 

7.319 

-.2046 

15 

1.4 

.0113 

.674 

-.0321 

.473 

-10. 3501 

-8.5878 

-15.2285 

-.07581 

.16976 

-10.555 

-.2049 

16 

1.5 

.0573 

.722 

.0345 

.545 

-22.6711 

-17.0053 

-6.4301 

-.16872 

.06866 

-22.775 

-.1039 

17 

1.6 

.0670 

.817 

.0491 

.658 

-23.8430 

-18.8982 

6.5720 

-.17860 

-.05885 

-23.785 

.0580 

18 

1.7 

.0243 

.884 

-.0217 

.753 

-13.9088 

-10.8324 

18.0802 

-.10522 

-.15221 

-13.7H 

.1978 

19 

1.8 

.0004 

.894 

-.0629 

.805 

1.8028 

3.2344 

26.3338 

.01196 

-.17032 

2.053 

.2502 

20 

1.9 

.0353 

.914 

.0020 

.871 

• 15.5842 

15.2862 

13.2823 

.11543 

-.10984 

15.764 

-1798 

21 

2.0 

.0706 

.994 

.0705 

.999 

21.1848 

21.0787 

.3580 

.15819 

-.00415 

21.207 

.0222 

22 

2.1 

.0474 

1.084 

.0225 

l.l4i 

16.6536 

17.1277 

-13.8916 

.12522 

.09490 

16.504 

-.1496 

23 

2.2 

.0044 

1.117 

-.0697 

1.229 

4.9112 

4.1234 

-21.3879 

.03794 

.14 266 

4.663 

-.2482 

24 

2-3 

.0135 

1.122 

-.0489 

1.291 

-8.0209 

-10.1541 

-18.2064 

-.05883 

.12177 

-8.250 

-.2291 

25' 

2.4 

.0594 

1.174 

.0591' 

1.413 

-16.1158 

-19.3557 

-7.3413 

-.11994 

.04793 

-16.219 

-.1032 

26 

2.5 

.0656 

1.269 

.0736 

1.591 

-16.1305 

-20.0720 

7.0345 

-.12083 

-.04042 

-16.060 

.0705 

27 

2.6 

.0216 

1.333 

-.0386 

1.734 

-8.9181 

-10.8053 

18.2719 

-.06694 

-.10187 

-8.637 

.2811 

28 

2.7, 

.0010 

1.342 

-.0928 

1.811 

1.7033 

3.3439 

19.9773 

.01174 

-.11129 

1.945 

.2417 

29 

2.8 

.0381 

1.365 

.0096 

1.913 

10.3942 

14.8375 

12.6147 

.07702 

-.06981 

10.555 

.1608 

30 

2.9 

.0710 

1.447 

.1030 

2.103 

43.4568 

19.4654 

.0845 

.10049 

-.00133 

13.468 

.0112 

31 

3.0 

.0447 

1.53 1 * 

.0250 

2.305 

10.1343 

14.9646 

-12.5013 

.07622 

.06024 

9.997 

-.1373 




E 

© 2 

= 4286.39139342 


E 0® = 

-35.532782573 





E 

®® 

= -228.77281991 


E 0® = 

51.02915144 





E 

©0 

= 35.627944979 


E 

© a = 

0.3725473211 





E 

©© 

= 0.754833829 


E © © ■ 

0.0542286301 





£ 

®@ 

= 1.43276343 


E © © = 

-0.014427701 





E 

® 2 

= 4l4o. 34148937 


E @ 2 

0.3777271561 





E 

®© 

= - 4.987891922 


E © © = 

-0.438887876 






XABiiE V 11 . — IRONY'S METHOD APPLIED TO EXAMPLE 
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O .2 .4 .6 .8 1.0 1.2 1.4 1.6 1.8 2.0 

Independent variable, t 

(b) Response q m (t) to the pulse F(t). 

Figure /.- Pulse input and associated response used in example I. 





Figure 2 - The variation with time of two quantities used in example 




Figure 3 Response to a step input used in example 



Response, qjt) Input, F(t) 



Independent variat 
(b) Response, qjt), to input, F(t), o i 




Figure 5- 
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Figure 6 -The variation with time of four quantities used in example JR. 





